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Abstract 

This paper discusses the response of a sequential and a Bayes 
minimum variance orbit determination system to errors in the force 
model. For the most part, the error is introduced into the central 
mass . In most instances the processors do not know of nor solve for the 
error. A half orbit of generated data at one minute intervals is used 
in all cases . 

It is shown that they both arrive at the same orbit even after 
the so-called "divergence" of the sequential processor. 

It is further shown that the residual patterns obtained in the 
final orbit contain frequencies which appear like the 6th order zonal 
harmonic even though the error was due to the zeroth order term. 




Finally, when this philosophy is applied in a manner similar tc 
the "raascon" reduction technique, a residual pattern is obtained which 
bears a striking resemblance to the lunar orbiter residual patterns 
although in fact no such "mascons" exist in the process. Hence, the 
technique of deriving "mascons’' from a differentiation of residual 


patterns is suspect . 



RESPONSE OF ORBIT DETERMINATION SYSTEMS TO MODEL ERRORS 


INTRODUCTION 

There are two fundamentally different orbit determination systems, 
namely, least squares or batch processing and sequential, usually 
Kalman processing. It is well known that under ideal conditions i.e. 
if all the assumptions made in their derivations are met; they converge 
to the same solution. In orbit determination, this means they find the 
same trajectory. 

Least squares has been in use for generations and is currently the 
most accepted for scientific analysis. Sequential processing has only 
recently been introduced and has seen only limited use on practical 
problems . 

The mathematical development of both methods has been widely pub- 
lished and wil] not therefore be repeated here. Occasionally the equa- 
tions will be used for the purpose of clarifying the discussion. 

Applying either of these techniques to orbit determination is a 
very complex task which, in general, tends to mask, if not obliterate 
the physical picture of what is taking place. Often they have been 
used without rigor i .e . in violation of one or more of the assumptions 
made in their derivation. It is no wonder that strange results occur. 
There exists a vast amount of literature on how these systems respond 
in the ideal, theoretical world but little exists which demonstrates 
how they respond to things approaching the real world. This is partly 
because in the real world the true or actual orbit is never known. At 
beet then, one can attempt to simulate certain features of the real 


1 


world and test the system responses under controlled conditions. This 
work represents at least an initial step in this direction. 

Two computer programs are used in this study, which are identical 
in every respect except in the statistical subroutine . One uses a Bayes 
batch processor, the other the Kalman sequential processor. In the 
latter, the capability exists to first scan the data to examine residual 
patterns before the differential correction pass and again to scan 
after the differential correction pass . 

In this study ideal observations were generated for a single half 
orbit of a point mass around the earth. The data are range and range- 
rate observations made from a single tracking station located at a 
fixed point in inertial space at about the 2k hour altitude. Data was 
generated at one minute intervals using the SAO '66 gravitational model. 
The orbit being tracKed had the following initial osculating elements: 

Epoch 

a * 2.099 ER 

e =* 0.500 

i = 100?00 

ou * 09003 

fi * 0?0 

t =» 1 .289 hrs . 

P 

The pericenter was an the equator and so was the tracking station. 
The noise added to the ideal data by a random number generator was 5*0 
meters in range and about 2.0 mm/sec in range-rate. 
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DIVERGENCE 


The first phenomenon to be discussed is that of the so-called 
’’Divergence" of the sequential processor. In Figure 1, the residuals, 
i.e. the observed minus computed observations, are shown for a sequen- 
tial orbit determination pass through the data. In this case, the 
central mass is in error by one in the fourth significant figure. This 
large error was selected so that the error was well above the noise . 
There has been no attempt to solve or account for this error in the 
statistical computations. A large and arbitrary a priori matrix has 
been used (Diagonal, cartesian, 1. Km each position coordinate and 
1. ra/s velocity) and with this low noise data, its influence is insigni- 
ficant after the first six observations. The upper curve shows the 
residual pattern when only range data are used and the obvious diver- 
gence occurs at about 0.9 hours. 

The lower curve shows the behavior when only range-rate data is 
used with divergence occurring at 1.2 hours. The time of divergence is 
strictly dependent on the noise level assigned to the data but the fact 
that they both occur at about 1. hour indicates that there is about the 
same information content in both observation lyp^s . 

In Figure 2, the same case is shown but both the observation types 
are used in the determination. Now divergence begins at 0.2 hours; of 
course the observation types are assumed independent. 

An examination of the orbits at the point where divergence occurs 
shows it to be the point where the covariance matrix can no longer 
represent the error, i.e. the model has broken dowh. This may not be a 
great surprise but continuing to process the data to 154 minutes did 
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DIVERGENCE 

SEQUENTIAL ORBIT DETERMINATION 
OBSERVED fi =19.9094165 
COMPUTED li - 19.9194165 
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SEQUENTIAL ORBIT DETERMINATION 
RAN6E AND RANGE-RATE OBSERVATIONS 
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no further damage to the orbit. The covariance matrix, however did get 
less and less realistic. All three of the cases obtained the same 
orbit to four significant figures which of course coincides witn the 
model error imposed. The point to emphasise, however, is that diver- 
gence tells you not only when your model has failed but in what 
significant figure the error is occurring. 

EQUIVALENCE 

The Bayes processor was applied to the same cases and in the next 
four figures the final scans of both processors are shown, i .e . there 
is no differential correction being made in these passes. Figure 3 is 
the case where only range data is used. Obviously, they htve both 
obtained the same orbits. In fact they agree to five or nr re signifi- 
cant figures. You may note some rather strange residual patterns. 

These will be discussed later. Figure 4 is the case where range-rate 
only was used. Figure 5 shows the range -residuals when both data types 
were used in the orbit determination and Figure 6 the range-rate. 

A point should be emphasized here that the Bayes does not show an 
obvious place where the model breaks down. ‘They both, however have 
matched covariance matrices after convergence has occurred. The Epoch 
Matrix of the Bayes can be propagated to t-raax to compare with the 
Kalman or the Kalman may be propagated back to Epoch witn the same 
agreement . 













BAYES-SEQUENTIAL EQUIVALENCE 

RANGE-RATF. OBSERVATIONS ONLY 
OBSERVED fL = 19.9094165 
COMPUTED 19.9194165 
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MINUTES AFTER EPOCH 





BAYES-SEQUENTIAL EQUIVALENCE 
RANGE AND RANGE-RATE OBSERVATIONS 
OBSERVED 11 = 19.9094165 
COMPUTED u = 19.9194165 
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RANGE AND RANGE-RATE OBSERVATIONS 
OBSERVED 19.9094165 
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Figure 6 
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MINUTES AFTER EPOCH 



MODEL ERROR EFFECT 


What happens to the system when the error in the model is accounted 
for? You may recall that in a sequential processor the state correction 
is given by Aa = L Ay where Acr is the correction to the state, Ay is the 
residual, i .e . observed minus computed observation and L is the Kalman 
gain matrix. The usual notation of Schmidt has been left to apply to 
the cartesian state. The notation here is one step up the alphabet for 
the alpha state hence the L rather than K for the Kalman gain matrix. 
When the state consists only of the six parameters 

L = or Y _1 

where N is the sensitivity of the observation to th? state dy/da, Q the 
covariance matrix and 

Y = NQN' + e " 2 

where e 2 is the white noise variance. 

/ Q c 

If now the state is augmented by we have ^ ^ ^ j for the 

augmented state covariance. By accounting for the error we mean com- 
puting the augmented covariance but not updating nor decrementing its 
covariance. This is the procedure given by Schmidt. The gain then 
becomes 

L = (QN' + C F ' ) Y ” 1 

Y - NQN' + pw F' + NC F' + FC ’N' + e " 2 

b pi jj, 

where F is the sensitivity of the observation to the new state variable 

& 


The next four figures show what happens to the sequential proces- 
sor during orbit deteriri nation when the bias is included. 

Figure 7 shows the range without the bias matrix in the upper 
curve when both observation types are used. The lower curve shows how 
it behaves when the bias is included. The divergence of course disap- 
pears and the orbit continues to follow the observations . 

Figure 8 is the range -rate . 

Figure 9 shows the patterns when only range data are used; 

Figure 10 when only range-rate data are used. 

The Bayes processor does not yet have this capability, but with 
the agreement we have already obtained, we have no reason to suspect 
that it will not get the same results as those ohown in the next four 
figures . 

Figure 11 shows the final range scan when both range and range-rate 
data were used. Obviously, the orbit fits the latter part of the data 
for the span of the model validity. 

Figure 12 is the range-rate data for the same case. 

Figure 13 shows what happens when range only is used and Figure 14 
when range-rate is used by itself. 

The advantage of using unsolved for biases should be apparent. 

The problem of modeling increases without limit as we proceed to higner 
and higher accuracies . Using unsolved for biases permits us to define 
a reasonable cutoff point at some finite accuracy and at some manageable 
computing complexity level. This point should vary depending on the 
purpose of the end product. Needless to say, prediction is always 
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MODEL ERROR EFFECT 

SEQUENTIAL ORBIT DETERMINATION 
RANGE AND RANGE-RATE OBSERVATIONS 
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MODEL ERROR EFFECT 
SEQUENTIAL ORBIT DETERMINATION 
RANGE OBSERVATIONS 
OBSERVED fi= 19.9094165 
COMPUTED M= 19.9194165 
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SEQUENTIAL ORBIT DETERMINATION 
RANGE-RATE OBSERVATIONS 
OBSERVED fL = 19 9094165 
COMPUTED 19.9194165 
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MINUTES AFTER EPOCH 




MODEL ERROR EFFECT 
SEQUENTIAL FINAL SCAN 
RANGE AND RANGE-RATE OBSERVATIONS 
OBSERVED u= 19.9094165 
COMPUTED/t =19.9194165 
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MODEL ERROR EFFECT 
SEQUENTIAL FINAL SCAN 
RANGE AND RANGE-RATE OBSERVATIONS 
OBSERVED u= 19.9094165 
COMPUTED^. = 19.9194165 
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MINUTES AFTER EPOCH 




SEQUENTIAL FINAL SCAN 
RANGE OBSERVATIONS 
OBSERVED li 19.9094165 
COMPUTED 19.9194165 
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Figure 14 
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model limited, therefore, prediction needs will dictate the required 
model accuracy. Beyond that the model should be cutoff with unsolved 
for biases. Only through this type of process can observations and 
orbit determination systems be utilized to their maximum efficiency. 


RESIDUAL FREQUENCIES 


The appearance of the final scan residuals led us to providing a 
preliminary scan capability in the program in order that we may study 
the complete evolution of the orbit determination process. 

Figure 15 shows the comparison of preliminary and final scan for 
the range only observation case. It is clear, that the orbit determina- 
tion process has introduced some spurious frequencies in the final 
residual pattern which could easily be mistaken for an error in J2,o or 
J 4 ,0 or sorr. . other coefficient rather than the central mass . We have 
suspected this for a long time but now there can be no doubt . In 
Figure 1 6 we see the same misfortune befalling the range-rate. It is 
at present only conjecture but we are quite confident that if we 
enlarge the state the situation will be far worse . 

There is some hope for this dilemma as we shall see in the next 
two figures when range and range-rate observations are used. Figure 17 
shows the range. Now the residual pattern is more like the error source 
that produced it, at least in the number of wave lengths and Figure 18 


agrees showing the range - rate . This is a strong case for at least two 
and possibly more independent data types . It also explains to sane 
extent why different data types find different orbits . Finally, it 
warns us to be cautious in the analysis of frequency content of residual 
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RESIDUAL FREQUENCIES 

SEQUENTIAL PROCESSING 
RANGE AND RANGE-RATE OBSERVATIONS 
OBSERVED a = 19.9094165 
COMPUTED u = 19.9194165 
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patterns when a single observation type is used or dominates the 
determination, i .e . far greater accuracy than any others. 

Having obtained these results and noting the lunar orbiter dif- 
ficulties we decided to simulate with the earth, a crude approximation 
to the raascon processing. 

In Figure 19, we see the residual pattern for range, when both 
range and range-rate were used as observations. In all of this series 
we used the same data as before but now, the error we have imposed is 
the use of a spherical earth for the orbit determination. 

In Figure 20 the range-rate residuals are shown — nothing alarming 
in either of these . 

In Figure 21 we see what the residuals are like when range only is 
used and now strange things are happening. 

In Figure 22 we see a very familiar pattern when range-rate is used 
alone. Before you jump to the conclusion that this is a fault of sequen- 
tial processing look at Figure 23 which compares sequential and Bayes 
with range-rate only. We should further say that all the runs, with the 
exception of the unsolved for bias have been verified using the Bayes. 

In the final Figure 24, we further demonstrate the comparison by 
showing the lunar orbiter residuals and ours on the same figure . 

I do not have the figure, but I can say that divergence occurred 
prior to reachixig perigee in the determination pass . 

We have of course only demonstrated a striking similarity on a 
single nass but we feel it certainly warrants further investigation. 

We offer no explanation for the high correlation of Mascons with surface 
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features we prefer to leave that to the Mascon investigators. 

In conclusion, we t>-ink we have begun a very enlightening study 
for those of us in the orbit determination business. This was begun 
as a prelude to the acquisition of continuous tracking of earth satel- 
lites from orbiting tracking stations. However, as complex as orbit 
determination models are, there is a large amount of room for other 



investigators to follow analogous research and we would like to encour- 
age as many of you as possible to do so. Our programs though not 
designed for export, are available in as is condition. 

Our sequential processor has an error analysis mode which we have 


been using to check out error analysis programs. Let it suffice to say, 


that there should be a requirement placed on all vendors of error 




analysis programs to include the generation of observations and dif- 


ferential correction in order to verify that it is properly including 
ie . ’feet of each error source. There is no better method of checking 
such a program. Test cases should include a separate one for each 
error source and all combinations and permutations thereof. 
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